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' Recent data on the electroproduction of hadrons on nuclei are discussed. The 

. effects of the nuclear absorption are investigated using the Lund model for electro- 
production on nucleons. A simple geometrical model with a minimal number of 

Q^' assumptions and free parameters is shown to describe the data reasonably well. 
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1 Introduction 



Electroproduction of hadrons on the nuclei is considered since long as a favourite 
process to investigate the nuclear absorption effects. On one side, the initial production 
process may be safely assumed to occur on a single nucleon and to be the same as in the 
electroproduction on a nucleon (proton or neutron) in the deuteron, since any screening 
corrections are negligible for deeply virtual photons exchanged in the process. On the 
other hand, there is a marked difference in the spectra of hadrons from these two processes. 
The obvious effect leading to such a difference is the secondary interaction of the produced 
hadron on its way out from the nucleus [1]. 

Recently, interesting data from the HERMES collaboration were published, enlarging 
our knowledge of the electroproduction on N, Kr and Xe nuclei both for the single spectra 
[2] and for the two hadron systems [3]. The direct comparison of data on hydrogen and 
on heavier nuclei allows to investigate the ratios of spectra for which many systematic 
uncertainties cancel. 

We present a comparison of these data with a simple model based on the PYTHIA 
code [4] for the electroproduction on nucleons and the geometrical scheme for calculating 
the absorption effects. Obviously, we are not the first to consider a model of this type. The 
unique feature of the Lund description of the production processes is the definite space- 
time picture associated with the generation of hadrons [5] . Therefore it was extensively 
used to describe the absorption effects in the hadroproduction on nuclei [6, 7, 8, 9, 10, 11]. 
Also recently various versions of the Lund models were exploited to analyse the HERMES 
data [12, 13, 14, 15]. However in these papers the geometrical picture of production from 
the one dimensional Lund model was carried directly into the 3D case, and/or coupled to 
the comphcated description of interactions with "the remaining part" of the nucleus (e.g. 
the system of transport equations). It is thus difficult to decide which piece of the model 
is responsible for the discrepancies with data. 

We discuss a very simple picture, in which only the obvious part of the Lund space- 
time development is used, and we supplement it with (equally obvious) pure absorptive 
effects. Clearly, such a model has a limited apphcability. However, we feel that the 
determination of the limits of applicability for a maximally simplified model is useful for 
any future extensions and more sophisticated descriptions. Moreover, it is instructive to 
check which features of the data can be described satisfactorily already in such a simple 
model and thus, which predictions seem to not depend on any details of the description. 
We restrict ourselves to the use of hadronic (and not partonic) degrees of freedom, since 
we discuss the low energy data for which the typical values are small. 

In the next section we present our model in detail. Then we discuss the results for the 
single particle spectra and the comparison with recent data on hadron pairs. We conclude 
with a short summary. 

2 The model 

We are using the Monte Carlo generator PYTHIA 6.203 [4] and generate more than 
a quarter milion of events per each nucleus, with a correct proportion of ep or en events 
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(the experimental data are shown as ratios of the distributions from heavier nuclei and 
from the deuteron). All the kinematical cuts from two HERMES sets of data are applied, 
either by setting the proper values of PYTHIA parameters, or explicitly in the program 
for the event analysis. 

We supplement the ordinary information provided by PYTHIA for each event by the 
values of one extra parameter from the generating algorithm. It is the GAM(3) parameter, 
set for each string break in the PYSTRF procedure, denoting the proper time tq (time 
measured in the string rest frame) between the string formation and its break. This time, 
corrected for the Lorentz dilatation, is used to calculate the distance between the string 
formation and string breaking point in the nucleus rest frame 

«/orm = 'TQVstrlstr- 

Then it is attributed as the lower limit for the formation length to the hadron originating 
from this break. We allot this length for the charged "stable" hadrons, i.e. pions, kaons 
and (anti)protons born directly from the string breaks as well as for the decay products 
("children" and "grandchildren") of the hadrons born directly. For the "last" hadron 
(or the pair of hadrons) of the event, which combines the "leftover" partons born in the 
former breaks, the last calculated proper time of the break is used.^ We checked that this 
prescription attributes the formation length to almost all the charged hadrons produced. 

The assumption that the string break point is the "birthplace" of a hadron is known 
to lead to an overestimation of the absorption effects. Thus the usual geometrical con- 
siderations of the string breaking and hadron formation in the Lund model [5] suggest 
that the formation length should include the distance between the string break and the 
"formation point", the first intersection of the lines of partons "born" in the neighbour- 
ing breaks [15]. Therefore one distinguishes often between two formation lengths: the 
"production length" and the "yo-yo length" [8]. The later choice was found repeatedly to 
yield too small absorption effects, unless significant absorption cross section for partons 
is assumed. Moreover, such an intersection occurs in fact only in the 1-1-1 dimensional 
Lund model; in (3 + 1)D the lines usually do not meet. 

Thus we keep the idea of a "two stage" formation of hadrons, but consider other ways 
of accounting for the second stage (rearrangement of partons to form a hadron). We tried 
two simple choices. In the first case a constant "recombination proper time" (of the 
order of 0.1 1 fm/c in the string rest frame) is assumed, leading to the extra piece of 
formation length proportional to the string velocity and the Lorentz factor 

Sform = (To + Tr)Vstr1str- 

In the second case for each hadron a "hadronization proper time" Th is assumed, and 
subsequently dilatated by a hadron Lorentz factor 

^form — 

We found that the last choice overestimates significantly the formation length, leading 
to the underestimation of the absorption effects for fast hadrons (the ratios of spectra 

^Thus typically each value of the formation length is attributed to more than one charged hadron. 
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do not decrease for high hadron momenta). Thus this case is not considered in the 
following. Note that in both cases we assume no absorption until the hadron is formed 
(the absorption cross section is zero for the state existing between the string break and 
the hadron formation). 

For the string creation point we generate the values of the distance from the centre of 
nucleus and the azimuthal angle 9 assuming a constant density in volume, i.e. choosing 
randomly from the flat distributions of Tq between and i?^, and cosO between —1 and 
1. Here R is the nuclear radius, calculated for each A from the standard formula (in fm 
units) 

R{A) = 1.12A^/3 _ 0.86^-^/1 

To estimate the absorption effects we calculated first the geometrical path from the 
creation point to the nuclear surface, i.e. 

Sgeom = tqcosO + J R"^ - r^sin^O. 




''form 



form 



geom 



Figure 1: The schematic picture of the string formation, string breaking and hadron formation 
inside a nucleus. The distances are measured from the string formation point. 

Then we evaluate the absorption length L for hadrons inside nucleus from the absorp- 
tion cross section a and nuclear density p, L = {ap)~^. The inelastic cross sections in the 
energy range corresponding to the HERMES conditions are about 50 mb for antiprotons, 
30 mb for protons, 20 mb for pions and negative kaons and about 12 mb for positive 
kaons. This yields the values of absorption length of 1.2 -^ 5 fm. With these values we 
may calculate the survival probability 




geom 



/ L) if Sg^^^ Sgeom 

if s' < 
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which is used as the weight to produce the spectra from nuclei. The ratio of spectra 
with weights and without weights (i.e. with weight equal one) is to be compared with 
the experimental ratios of spectra from nuclei and the properly normalized spectra from 
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nucleons. The geometrical picture showing the quantities defined above is presented in 
Fig.l for the case of a hadron formed inside the nucleus. 

In these calculations the nuclei were modelled by a very simple description, the so- 
called "hard sphere", in which the nuclear density is constant within a sphere of radius 
R, and zero outside. In a more realistic model the nuclear density depends on the radial 
coordinate e.g. by the standard (" Saxon- Woods" ) formula 

p{r) 



l + exp[{r- R)/Ay 

where po — 0.17 fm~^ and A = 0.54 fm. To check if the geometrical approximation is 
reliable we have repeated the calculations with s'g^g^ replaced by s^^^ 



4// = / 



oo 



d 



1 -|- exp 



A 

where d = —TqCosO + Sform- Another modification was the selection of the string creation 
point not just from the inside of the sphere, but with tq weighted with the Saxon- Woods 
distribution. 

The ratios of ^-distributions resulting from the weights calculated in the discussed 
cases for all nuclei are hardly distinguishable. Since the generation of values from a 
non-uniform distribution is more time consuming, we use in the following the model with 
variable nuclear density, but with the simple generation of the Tq values. 



3 Absorption for the single hadron spectra 

The HERMES collaboration has published precise data for single charged hadron 
electroproduction on the dcuteron, nitrogen and krypton gas targets for the positron 
energy of 27.6 GeV [2]. The relative energy in the lab frame z = Eh/u was used to define 
the spectra in the range 0.02 < z < 1.0. The ratio of the properly normalized spectra 

refiects the nuclear effects for the light (nitrogen, A = 14) and medium-size (krypton, A = 
84) nuclei. The kinematical cuts for the momentum transfer squared, > 1.0 GeV/c^, 
hadronic invariant mass W > 2 GeV and the positron relative energy loss y = v / E < 0.85 
were used. 

These results were compared in the HERMES paper [2] with two absorption models 
[12, 13]. One of them [12] visibly overestimates the absorption effects for nitrogen, and 
the other one [13] predicts a too steep 2;-dependence for the krypton data. Obviously, 
none of the models describes the excess of hadrons at small z, which is probably due to 
the production in secondary collisions. 

Wc compare the same data with the predictions of our model specified in the previous 
section. For the formation length defined by the breaking point of the string, Sform = 
TojstrVstr, the absorption effects are too strong, especially for the krypton data. Thus it 
seems necessary to introduce the "recombination proper time", as described above. 
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Figure 2: The experimental ratio of the charged hadrons z-spectra for nitrogen and deuterium 
[2] (stars) compared with the model calculations for = 0.7 fm/c (broken line) and for = 0.8 
fm/c (solid line). 



The value of was not fitted to any particular set of data, but we found that the 
value of Tr = 0.7 0.8 fm/c gives a good description of the data for nitrogen for z > 0.1 
(Fig. 2). The dependence of the results on this value is not very strong. In the following 
we use only r,. = 0.8 fm, which describes better the data for krypton for z > 0.3 (Fig.3). 

The good agreement with data of a simple model with only one free parameter for 
^ > 0.3 is remarkable. One should note that both the z-dependence and A-dependence of 
data seems to be correctly described using the geometrical picture provided by the Lund 
model. Obviously, it is desirable to get similar data for more than two values of A, e.g. 
for xenon, for which the predictions are also shown in Fig.3. 

It is obvious that we cannot expect an agreement of the model calculations with the 
data for small values of z. We have considered merely the absorption effects (and get in 
fact the results quite similar to other authors cited in the HERMES paper [2], although 
our model gives a better description of the data for large values of z). We expect that 
the secondary collisions of the hadrons originating from the string breaking with other 
nucleons produce extra hadrons. This enhances the spectra for low z values and may even 
result in the values of Rh above one, as seen in the first bin of the nitrogen data. Such an 
effect may be directly incorporated into the model [16], or combined with the absorption 
effects in the transport equations [14, 15]. However, such refinements are rather arbitrary 
and involve many extra parameters. Therefore we feel that it is useful to establish the 
range of z for which a simple absorption model does not require any such modifications. 
It seems to be almost the full z range for light nuclei, as nitrogen (z > 0.1), and visibly 
less for heavier nuclei, as krypton (z > 0.3). It would be interesting to check if the same 
range applies as well for other nuclei, e.g. xenon. 

Till now we have discussed the data for all charged hadrons. The HERMES experiment 
is able to identify particles and the corresponding data were also published [2]. The 
data for pions resemble closely those for all charged hadrons, as pions are most copiously 
produced. The data for negative kaons do not agree with the predictions of the absorption 
model [12]. In particular, this model predicts almost negligible difference between the 
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Figure 3: The experimental ratio of the charged hadrons z-spectra for krypton and deuterium 
[2] (stars) compared with the model calculations for Ty = 0.8 fm/c (solid line). The model 
calculations are shown also for xenon (broken line). 

positive and negative kaons, in clear disagreement with data. 

In our model one may expect different absorption effects for pions and kaons for two 
reasons. First, the absorption cross sections are different, as discused in the previous sec- 
tion. Second, the different masses may lead to different values of the only free parameter, 
the recombination time r^. We found that a reasonable agreement with data for kaons 
at z > 0.35 is achieved for Tr = 0.4 fm/c, half of the value for pions. Let us note that 
by the uncertainty relation these values correspond to the energies of about 0.3 GeV for 
pions and 0.6 GeV for kaons, close to their typical transverse masses. The data are shown 
together with our predictions in Fig.4. 
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Figure 4: The ratio of normalized yields of K^-s (solid squares) and K~-s (open squares) for 
Kr and D nuclei [2]. Our model predictions are shown as broken and solid lines, respectively. 
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4 Two hadron spectra 



Recently, the HERMES collaboration has extended their previous data by publishing 
precise results on the conditional z-spectra for charged hadron electroproduction on the 
deuteron, nitrogen, krypton and xenon gas targets for the same positron energy of 27.6 
GeV [3]. Kinematical cuts were the same as for the single spectra [2]. The events were 
selected by the condition that one of the charged hadrons has zi > 0.5. The ratio of the 
properly normalized spectra for other hadrons 

^2/^1^2) - [ J^^^ )J[ J^^^ 

(of the same charge sign as the selected one) were shown for 0.05 < Z2 < 0.5. 

The results were compared with the predictions of absorption models [12, 13] and 
with the predictions of a model including the production processes by means of the BUU 
transport equations [14]. The absorption models overestimate the absorption effects, 
especially for heavier nuclei, in all the considered 2;-range. The transport model fares 
better, but it does not reflect two most striking features of the data: the similarity of 
the results for all the nuclei and the apparent disappearance of the absorption effects for 
Z2 0.5. 
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Figure 5: The experimental ratio of the conditional charged hadrons z-spectra for nitrogen and 
deuterium [3] (stars with errorbars) compared with the model calculations for Tj. = 0.8 fm/c 
(broken line) and the results of the event generation without proper correlations (solid line). 

We compare the data with the same absorption model which we used to describe 
the single spectra (with the value of tr = 0.8 fm). Let us note that a simple selection of 
PYTHIA events for ep and en collisions with the value of Zi > 0.5 for one of the hadrons is 
not sufficient to model the nuclear collisions. Due to the strong nuclear absorption of the 
fast "trigger hadron" the cut on z discriminates against the selection of a hit nucleon in 
the "front part" of the nucleus. Therefore, on average, the other hadrons will be absorbed 
less strongly than in a typical event. 

In the Lund models this effect is obviously present. The observation of a fast hadron in 
electroproduction on heavy nucleus enhances the probability that the absorption of all the 
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hadrons in this event is below the average. This happens because the absorption depends 
both on Sform (which may be different for each hadron) and Sgeom, which is common for 
all the hadrons from a given event. Therefore a proper procedure to generate events for 
the conditional ^-distributions is to calculate for each event with a fast hadron {zi > 0.5) 
the probability P„a that this hadron is not absorbed before leaving the nucleus. Then one 
generates a random number r from the range -h 1 and accepts the event only if r < P„a. 
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Figure 6: The experimental ratio of the conditional charged hadrons z-spectra for krypton and 
deuterium [3] (stars with errorbars) compared with the model calculations for Tr = 0.8 fm/c 
(broken line) and the results of the event generation without proper correlations (solid line). 
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Figure 7: The experimental ratio of the conditional charged hadrons z-spectra for xenon and 
deuterium [3] (stars with errorbars) compared with the model calculations for Tr = 0.8 fm/c 
(broken line) and the results of the event generation without proper correlations (solid line). 

In figures 5-7 we show both the predictions of our model with the correlation taken into 
account, as described above, and without this correlation (showing just the absorption 
for the hadrons other than the fastest one). We see that the nitrogen data show stronger 
absorption than expected (and stronger than seen in the single spectra, c.f. Fig. 2), but 
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the disagreement is in the worst case on the level of two standard deviations. The data for 
krypton and xenon are marginally consistent with the model for Z2 > 0.2 (Kr) or Z2 > 0.3 
(Xe). The important point is that the model with correlations is much closer to the data. 

Let us remind once more that for lower z the production processes (neglected in 
our model) are probably important. The results for the model neglecting correlations 
overestimate significantly the absorption effects for heavier nuclei. They are quite similar 
to the results of the absorption models shown together with the HERMES data [3] . 

5 Conclusions 

We have investigated the electroproduction of hadrons inside the nuclei using the 
PYTHIA generator. The results from the recent HERMES experiment [2, 3] are compared 
with the calculations done using a very simple scheme for the absorption effects. The 
ratios of single particle spectra for all charged hadrons are quite well described for the 
values of z above 0.1 (N) or 0.3 (Kr). For the lower z values our model predicts a 
flat 2;-dependence (as any absorption model), in disagreement with data. The flavour 
dependence of the absorption effects in single spectra is resonably well described assuming 
that the recombination time is proportional to the inverse of the hadron transverse mass. 
The ratios of the conditional spectra for the events containing one hadron with large 
z are also consistent with data for the values of z above 0.1 (N) or 0.3 (Kr, Xe). In 
particular, the similarity of the conditional spectra for all nuclei and the increase of ratios 
for maximal possible Z2 seems to be well understood in the absorption model with the 
PYTHIA generation mechanism. These effects result from the correlation induced by the 
bias in the position of a string creation for the events with a very fast hadron. 

Our results suggest that the space-time production picture in the Lund model coupled 
with a simple absorption model leads to a satisfactory description of the large part of the 
experimental data. 
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